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We generalize the Wannier interpolation of the electron-phonon matrix elements to the case of 
polar-optical coupling in polar semiconductors. We verify our methodological developments against 
experiments, by calculating the widths of the electronic bands due to electron-phonon scattering 
in GaAs, the prototype polar semiconductor. The calculated widths are then used to estimate the 
broadenings of excitons at critical points in GaAs and the electron-phonon relaxation times of hot 
electrons. Our findings are in good agreement with available experimental data. Finally, we demon¬ 
strate that while the Frohlich interaction is the dominant scattering process for electrons/holes close 
to the valley minima, in agreement with low-field transport results, at higher energies, the intervalley 
scattering dominates the relaxation dynamics of hot electrons or holes. The capability of interpo¬ 
lating the polar-optical coupling opens new perspectives in the calculation of optical absorption and 
transport properties in semiconductors and thermoelectrics. 
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I. INTRODUCTION 

Electron-phonon coupling plays a fundamental role 
in the relaxation of photoexcited electrons, thus af¬ 
fecting the performance of photovoltaic! and other 
semiconductor-based device!! In many cases, the 
electron-phonon coupling determines the magnitude of 
the lifetimes of electronic states inside the band gap, and 
the widths of the corresponding absorption peaks!!. It 
is also responsible for shifts and broadenings of the inter¬ 
band critical points in semiconductors!! The interpre¬ 
tation and analysis of the (ultra)fast dynamics of relax¬ 
ation of photoexcited electrons is particularly difficult, 
because many relaxation processes are present simulta¬ 
neously, and disentangling their respective contributions 
requires ad hoc information on their relative importance 
and order of magnitude! 

At the same time, ab initio calculations provide 
an effective tool to estimate electron-phonon coupling 
strength in metals! — and semimetals like graphene!!— 
or bismut h 17 : 18 . In the case of semiconductors, the pre¬ 
dictive capability of calculations based on density func¬ 
tional perturbation theory (DFPT) 19 : 20 for the electron- 
phonon matrix elements has been demonstrated in a 
number of semiconductors!!—, alloy s 24 : 25 and nanos¬ 
tructure s !—. 

Recently, a method to interpolate the electron-phonon 
coupling matrix elements using Wannier fuctions has 
been introduce d 11 : 27 : 28 , providing a computationally ef¬ 
ficient method to calculate electron-phonon matrix ele¬ 
ments on extremely fine grids in the Brillouin zone (BZ) 
of metals. This has proved to be crucial to predict various 
material properties such as for example nonadiabaticity! 
or superconductivit} : 29 ' 30 . 


The metho d 11 : 27 has also been used to increase the pre¬ 
cision of integrals related to electron-phonon scattering 
times and scattering rates in semiconductors!, being, 
however, limited to nonpolar semiconductors. Indeed, 
in polar semiconductors, a long-wavelength longitudinal 
optical (LO) phonon induces an electric field, and the in¬ 
teraction of electrons with this macroscopic electric field 
- the polar-optical coupling or Frohlich interaction - is 
divergent when the phonon wave vector tends to zero. 
As the electron-phonon matrix elements related to the 
Frohlich interaction are not localized in the Wannier ba¬ 
sis, these matrix elements cannot be properly interpo¬ 
lated with the method of Refs. Em 

In this work, we extend the metho d! 1 :- 7 to take into 
account polar-optical coupling, and apply it to GaAs 
which is an archetype of a polar semiconductor. First, 
we present the theoretical background of our method, 
which we validate by comparing the Wannier interpolated 
electron-phonon matrix elements with that obtained by 
direct calculation within DFPT. Next, we calculate band 
broadenings due to the electron-phonon interaction for 
the highest valence and lowest conduction states. The 
calculated broadenings represent the total probability of 
the momentum relaxation of the hot electrons or holes 
due to electron-phonon interaction, in good agreement 
with recent pump-probe experiments— We analyze the 
role of the Frohlich interaction in the relaxation of excited 
electrons in GaAs, and we find that the Frohlich inter¬ 
action is responsible of the quasi-totality of the electron- 
phonon relaxation rates at low excitation energies, while 
representing only 10% of the electron-phonon relaxation 
rates for hot electrons or holes. The calculated data have 
then been used to estimate the broadenings of the E\ and 
E 2 critical points (CP) in GaAs. For E\, the calculated 
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broadening is in very good agreement with the experi¬ 
mental results of work [7! and with previous calculations^. 
Finally, for £ 2 , the calculated broadening are in satis¬ 
factory agreement with experimental results, in contrast 
with the previous calculation with the empirical pseu¬ 
dopotential method^. 


II. THEORY 

A. Electron-phonon matrix element 

The matrix element d^„(k, k + q) of the periodic part 
of the static and self-consistent response potential for 
a monochromatic perturbation of q wave vector, u qs , 
reads: 

d mn( k ^ k + q) = (H ^ SCF |k + qm), (1) 

where |kn) stands for the periodic part of the Bloch 
wave function of the initial electronic state, i.e. IV'k,™} = 
e lkr |kn)/v^Vfc- The vector k is the electronic wave vec¬ 
tor, q is the phonon wave vector, and n and m are the 
band numbers of the initial and final states. is the 
number of points in the k-grid on which V'k.n are gener¬ 
ated, the periodic part of the wave-function being nor¬ 
malized in the unit cell. u qs is the Fourier transform 
of the phonon displacement of atom s. The quantity 
SvscF/Suq S is the periodic part of the (static and self- 
consistent) response potential. 

The electron-phonon matrix element gf lm reads: 

<m( k - k +q) = 5I eS1/ ( q ) d ™"( k ’ k+q )/'v/ 2MsW q ,y - 

s 

We have used e sy (q) for the phonon eigenvector (s 
labels the atoms in the unit cell, v labels the phonon 
mode), uiqy is the phonon frequency, and M s is the atomic 
mass. 

In analogy with our previous works [33lf34l and the orig¬ 
inal workl35l the deformation potential for an individual 
transition is defined as a quantity proportional to the 
absolute value of the electron-phonon matrix element of 
eq. ©: 

J L( kk + q) = ^ 2P ^ UJ — \9nm(Kk + Cl)\. (3) 

Here, p is the mass density of the crystal, and H is the 
crystal volume. In the case of several initial and/or final 
electronic bands, we define the total deformation poten¬ 
tial as: 


Dfot = y/X nm (DZ m )*. (4) 

The Wannier interpolation of electron-phonon matrix 
elements was first introduced in Ref. GS Implementation 
of the Wannier interpolation procedure into the Quan¬ 
tum ESPRESSO package^, which we have used in this 


work, was described in Ref. [27l . In this work, we repeat 
only the part of the comprehensive description of Ref. 
that is necessary for the understanding of the extension 
to polar-optical coupling introduced in the next subsec¬ 
tion. 

A set of Wannier functions centered on site R are de¬ 
fined by the relation: 

|Rm) = ~4= y^e~ tkR I/ nrra (k)|V> kn ). (5) 

v N * k „ 

A transformation matrix, U mn { k), is determined by the 
Wannierization procedure (see Ref. [13). 

The matrix elements d(/ ra (k, k + q ) a re calculated 
within DFPT. As emphasized in Ref. [27), the periodic 
parts |kn) and |k+q?n) have to be exactly the same wave- 
functions used for the Wannierization procedure. This al¬ 
lows to fix their arbitrary phases appearing in the |kn)’s 
because of the numerical routine used for the diagonal- 
ization of matrixes containing complex numbers, or other 
numerical reasons^!. 

The matrix element d in the Wannier function basis is 
obtained by Fourier transform as: 

1 N k 

d(U(R,R L ) = ^EE e ' !kR+!qRL d^,„,(k + q,k) 

k.q m',n' 

( 6 ) 

with 


,(k + q,k) = b^ m ,(k + q)d^,„,(k + q,k)LT n ,„(k) 


( 7 ) 


Finally, when the localization conditions are verified 
on d^„(R, R L ) (see Ref. H3), one can obtain, by a slow 
Fourier transform, (k + q, k) with k and q being any 
points in Brillouin zone: 


d mn( k ~bq,k) (Nh'j 2 


EEE 




L R, m'n' 

Um'mOi + q)d*„/„/ (R, H L )U* n , (k) (8) 


In this work, GaAs is described within the local density 
approximation (LDA), and with the same pseudopoten¬ 
tials as in our previous works mui For the electronic 
density calculation, we used an energy cutoff value of 45 
Ry and a Monkhorst-Pack grid of 12 x 12 x 12 points in 
the BZ. The Wannier interpolation of the structure was 
carried out using 6x6x6 and 8x8x8 k-point grids cen¬ 
tered at r, with ten Wannier functions and 45 DFT Bloch 
wave functions. The large number of Wannier functions 
and DFT bands, as well as the rather dense k-point grids, 
are related to the costly disentanglement procedure nec¬ 
essary to satisfactorily reproduce the lowest conduction 
bands of GaAs. 

The same 6x6x6 and 8x8x8 grids centered at 
r were used as initial grids to calculate electron-phonon 
matrix elements within DFPT, which were then Wannier- 
interpolated using the interpolation method extended to 
polar-optical coupling described in the next paragraph. 
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FIG. 1: GaAs. Deformation potentials for the total contribu¬ 
tion of the three highest valence bands, see Eq. 0. The initial 
electronic state is at k = T, and the phonon wave vector q is 
varying along high symmetry lines in the BZ. Black line: di¬ 
rect DFPT calculation. Grey squares: Wannier-interpolation 
extended to polar-optical coupling. Circles: standard Wan- 
nier interpolation for metals and non-polar semiconductors. 


B. The long-range Frohlich interaction 


The Frohlich interaction is long-range, and thus the 
electron-phonon matrix elements for long-wavelength 
longitudinal optical phonons in polar materials are not 
localized in the real-space Wannier basis. A proof is given 
in Fig. [I] where the interpolation method of Refs. Illll27l 
is shown to fail for LO phonon as q —> 0. Indeed, as 
one can see from Fig. [1] at large |q| vectors, the charac¬ 
ter of the electron-phonon matrix elements is completely 
short-range for all phonon branches including the LO 
branch and is well reproduced by the standard interpo¬ 
lation method^. At small |q| vectors however, instead 
of the characteristic l/|q| behaviour, the interpolation 
method of Ref. [27] yields the same values of deformation 
potentials for the LO and TO branches. This is because, 
for the highest valence bands, the electron-phonon inter¬ 
action with long-wavelength LO phonons at k = T con¬ 
tains both long-range (Frohlich) and short-range contri¬ 
butions^. At vanishing |q|, the short-range contribution 
is the same for the LO and TO phonons. 

The absence of LO/TO splitting within a standard 
Wannier interpolation procedure is in analogy with the 
LO/TO splitting of phonon frequencies in polar semicon¬ 
ductors. Indeed, if dynamical matrices are interpolated 
with a real space cutof f 20 ' 38 , the LO/TO splitting of the 
phonon modes is absent, and the long-range part of the 
dynamical matrix needs to be subtracted before Fourier- 
interpolation into the real space and re-introduced after¬ 
wards in order to reproduce the LO/TO splittin g 20 ' 38 . 
We apply similar scheme in the case of electron-phonon 
matrix elements, with the nonlocal part of the electron- 
phonon interaction represented by the Vogl model de¬ 
scribed in the next paragraph. 



FIG. 2: GaAs. LO deformation potentials for the initial elec¬ 
tronic state with k = X and the phonon wave vector chang¬ 
ing along the (100) direction: q = (q x , 0, 0)2n/a, for q x > 
0.01. Panel A (top): deformation potentials calculated within 
DFPT for the lowest conduction band (dashed line) and two 
lowest valence bands (dot-dashed line, visible in panel B), 
compared with the deformation potentials calculated from the 
model of eq. 0 (solid line). Panel B (bottom): The difference 
between the deformation potentials calculated with DFPT, 
and the model of eq. 


C. The Vogl model 

The electron-phonon matrix element of the interac¬ 
tion of electrons with the macroscopic electric field in¬ 
duced by long-wavelength longitudinal optical phonons 
(the Frohlich interaction) was derived by Vogl in Ref. |H, 
for long-wavelength phonons (q —> 0). The leading con¬ 
tribution, in terms of ascending powers of q , to the in- 
traband electron-phonon matrix element, is given by the 
interaction with a dipole potential (oc jij) screened by 
the high-frequency dielectric tensor e*, (see Eq. 3.12 of 
Ref. 1391 1 . and reads: 

9diel(0.) ~ / ' -EE^- er(q)/v/2 M s uv 

q • £ °0 • q s y 

(9) 

The corresponding deformation potential then becomes: 

q) = ^^l!i,(q)l (10) 

In eq. (0, e is electronic charge, Zy Xs is the Born ef¬ 
fective charges tensor for atom s, and A, A' denotes the 
cartesian components. Among other approximations, it 
was assumed, in Ref. [H, that: 

(k + qm\kn) =S mn + 0{q 2 ). (11) 












Note that such a relation assumes a smooth and analytic 
relative phase relation among the |k, n) and |k + q, m) 
states. Such a requirement is satisfied by the phase choice 
giving the localised Wannier functions, but not by the 
arbitrary phase given by the diagonalisation procedure. 
Furthermore, eq. mu is the reason why the expres¬ 
sion m does not depend on the electronic wave vector k 
nor on the band indexes of the initial and final electronic 
states (which are assumed to be the same). Nevertheless, 
expression © describes well the asymptotic behaviour 
of the electron-phonon matrix elements in polar semi¬ 
conductors as q —l 0, as one can see in Fig. [2j where the 
behaviour of the electron-phonon matrix elements for the 
lowest conduction band (dashed line) and for the high¬ 
est valence bands (dot-dashed line) of GaAs calculated 
within DFPT are compared with the equations ([2]) and 
([Toll of the model for the Frohlich interaction along the 
(100) direction in the Brillouin zone (solid line). At large 
q x , the short-range character of the electron-phonon ma¬ 
trix elements is different for the conduction band and the 
valence bands, and cannot be described with the model 
of Eq. ©. At small q x , on the contrary, the asymptotic 
behaviour of the electron-phonon matrix elements for the 
valence bands and the conduction band becomes similar, 
and this behaviour is extremely well described by the 
model of eq. ©. 


D. Wannier interpolation extended to polar-optical 
coupling 

The method we propose in order to extend the inter¬ 
polation method of the electron-phonon matrix elements 
to polar semiconductors is similar to the one described 
in Ref. [38j for the interpolation of the force constants in 
polar materials. The idea is that the long-range con¬ 
tribution to the electron-phonon matrix elements, de¬ 
scribed with the model of eq. [9] has to be subtracted from 
the electron-phonon matrix elements before the Fourier 
transform to real space in eq. © and restored after the 
Fourier transform back to the reciprocal space in eq. ©. 

We use the Ewald sum in order to take into account 
the periodicity properties of the crystal and define: 

dt el ( q) = 

e -(q+G) 2 / 4 a 

(q + G) + G) Z «’ fe ' + G »'> 

A ' Jr 

( 12 ) 

Here, a is a convergence parameter (we used a = 5 (^y) 2 
in this work ), and the G are the reciprocal lattice vec¬ 
tors. The dielectric term, eq. urn depends only on the 
phonon wave vector q, and on material characteristics 
such as the Born effective charges Zy Xs and the dielectric 
constant eoo, which are calculated within linear response 

theory2£. 

Then we subtract d* ei (q) from the intraband Fourier 
transform of the deformation potential in the optimally 



FIG. 3: GaAs. Deformation potentials for the total contribu¬ 
tion of the three highest valence bands, see eq. multiplied 
by the modulus of the phonon wave vector |q|. The initial 
electronic state is at k = T, and the phonon wave vector q 
is varying along high symmetry lines in the BZ. Black line: 
direct DFPT calculation. Grey squares: our results obtained 
with the Wannier interpolation extended to the polar-optical 
coupling. 


smooth subspace, namely we define: 

Vm'n' ( k + q- k ) = dm'n' ( k + k ) “ <W ,n' df° l (q) 

(13) 


We then carry out the transformation in Eq. [S] on the 
matrix T>^, ra ,(k + q, k), 




EEE 


^kR+zqR,^ 


u„ 


L R, m'n' 

(k + q)^,„ ; (R,R L )C/„* n ,(k) (14) 


where I>^„,(R, Rl) is obtained via Eq. [Hlwith d s m , n , (k+ 
q, k) replaced by V s m , n ,{ k + q, k). 

Finally we add back d* e, (q) where now q is any 
phonon wave vector in the Brillouin zone, namely 


d^n( k + q- k ) = ^mn( k + q< k) + dff (q)£ m ,„ (15) 


In figure [T] the results obtained with the method of 
Wannier interpolation extended to polar-optical coupling 
are shown in grey squares. The behaviour of deformation 
potentials corresponding to LO phonons is well repro¬ 
duced by the Wannier interpolation extended to polar- 
optical coupling, in contrast with the standard Wannier 
interpolation method. 

In figure O we show the total deformation potentials 
for the three highest valence bands of GaAs (n, m = 
2,3,4), multiplied by the modulus of the phonon wave 
vector |q|, for all six phonon modes of GaAs. The crystal 
momentum of the initial electronic state was taken to be 
k = F, while the wave vector of the final electronic state 
k + q changes as the phonon wavector q varies along high 
symmetry lines in the BZ. We chose to multiply defor¬ 
mation potentials by the modulus of q, as, due to the 









Frohlich interaction, the deformation potential for the 
LO phonon tends to infinity as -j-ij and thus the values 
are very high close to T. In black are represented ref¬ 
erence DFPT values, and in grey squares are shown the 
deformation potentials which were interpolated using our 
Wannier interpolation method extended to polar-optical 
coupling. As one can see, the agreement between DFPT 
calculations and the Wannier-interpolated deformation 
potentials is excellent. The non-zero value of the defor¬ 
mation potential multiplied by q at q = F is due to the 
diverging LO-phonon Frohlich interaction, which is now 
properly described. 

In conclusion, the method of Wannier interpolation 
extended to polar-optical coupling yields interpolated 
electron-phonon matrix elements with the same preci¬ 
sion as the ’’standard” one at large phonon q vectors, 
but, in contrast to the standard method, it enables us 
to correctly describe the diverging LO-phonon Frohlich 
interaction at vanishing q. 

III. RESULTS 

A. Scattering rates and role of the Frohlich 
interaction 

The method of interpolation of the electron-phonon 
matrix elements in the Wannier space is necessary when 
one has to calculate integrals involving many q points, as 
it allows to significantly reduce the computational cost, 
compared to direct DFPT calculation. We have applied 
the method described in previous section to calculate the 
total probabilities of the electron-phonon scattering for 
an electron initially in the lowest conduction band of 
GaAs, and for a hole initially in the highest valence band 
of GaAs. 

In figure [J we show the full width at half maximum 
r due to the electron-phonon coupling, which was cal¬ 
culated for the lowest conduction band and the highest 
valence band of GaAs as a function of the k vector of the 
initial electronic state, at a temperature of 300 K: 

r nk = Y I d( l\9nn' (k, k ± q)| 2 

, JBZ 



Here, e n k are the electronic eigenenergies, and 7V q is 
the phonon occupation number which is described by the 
Bose-Einstein distribution function. Upper and lower 
symbols refer to absorbtion and emission, respectively. 

In practice, the delta function in eq. m was replaced 
by a Gaussian function in order to calculate numerically 
the integral in eq. m- The integration was performed 
on a 48 x 48 x 48 q-point grid in the BZ. The calculation 
was converged with respect to the Gaussian broadening 
starting from 15 meV broadening. As the pseudopoten¬ 
tials used here reproduce well the respective positions of 



FIG. 4: GaAs. Panel A: highest valence and lowest con¬ 
duction bands broadened by the electron-phonon interaction. 
The broadening is due to all electron-phonon interaction pro¬ 
cesses. For the sake of visibility, the calculated broadening Fk 
of the electronic bands has been multiplied by 5 in this panel. 
For panels B-D, the black solid line describes the behaviour 
of the initial electronic state in the highest valence band and 
the grey solid line is that of the initial state in the lowest 
conduction band. Panel B: broadening Fk due to electron- 
phonon interaction, as a function of the k vector of the initial 
electronic state. Panel C: electronic density of the final states 
allowed by conservation laws, as a function of the k vector 
of the initial electronic state. Panel D: relative contribution 
of the Frohlich interaction only to the total broadening. The 
calculations were done at T =300 K. 


the minima of the conduction band of GaAs^i, the Kohn- 
Sham band structure values at equillibrium were used to 
calculate the integral (fl6l) . 

As one can see, the broadening due to electron-phonon 
coupling varies from a few meV at the bottom of the con¬ 
duction band or at the top of the valence band (i.e. for k 
close to r point of the BZ), to several tenths of meV at 
high initial electron energies for the conduction states, 
or low initial hole energies for the valence band. The 
behaviour of the total electron-phonon scattering proba¬ 
bility is similar to the one of the density of the final elec¬ 
tronic states allowed by the energy and momentum con¬ 
servation laws (panel C), as the probability grows when 
more final states are available for the electron-phonon 
scattering. It is, however, not exactly the same, as the 
electron-phonon matrix elements are not constant over 












the Brillouin zone. 

The contribution due to the Frohlich interaction is of 
a few meV and does not change much over the Brillouin 
zone. It is, however, the dominant scattering process for 
the electron close to the bottom of the T or L valleys, 
and for the hole close to the top of the valence band, as 
one can see from the Panel D of Fig [4] Away from the 
band extrema, the intervalley electron-phonon scatter¬ 
ing mechanism rapidly becomes the dominant scattering 
mechanism. This result is in agreement with available lit¬ 
erature for low-field transport^. Indeed, at ambient tem¬ 
peratures, the Frohlich interaction is the dominant scat¬ 
tering mechanism which determines the low-field trans¬ 
port in GaAs^S. At high fields, however, the Frohlich 
interaction no longer plays the main role, and the inter¬ 
valley scattering is expected to determine the relaxation 
dynamics of electrons and/or holes, as can be deduced 
from Fig. [4j In this respect, GaAs behaviour in simi¬ 
lar to that of non-polar semiconductors, i.e. silicon or 
germanium^. 


B. Relaxation times related to electron-phonon 
coupling 

The widths of the electronic levels due to electron- 
phonon coupling presented on Fig. [I] can be used to es¬ 
timate the relaxation times of hot electrons related to 
electron-phonon scattering. Recently, relaxation time of 
hot electrons excited in the CB of GaAs close to T at 
excess energy e ex = 0.78 eV with respect to the CB bot¬ 
tom was found to be 22±3 fs at 293 K—. We find the 
electron-phonon scattering time Tk = to be 30 fs at 
e ex = 0.78 eV, in satisfactory agreement with the exper¬ 
imental result of Ref. Irj . 


C. Broadenings of some critical points 

The widths of the conduction and valence bands pre¬ 
sented in Fig. []] can be also used to estimate the 
temperature-dependent broadenings of excitons at criti¬ 
cal points in GaAs, attributed mostly to electron-phonon 
scattering. In principle, one should rely on the many- 
body excitonic wavefunction to obtain the excitonic life¬ 
time. The latter consists in a combination of prod¬ 
ucts of electron and hole quasiparticle (QP) wavefunc- 
tions linearly mixed through the exchange opera¬ 

tor - the electron-hole interaction 4 ^. Obtaining a precise 
knowledge of the QP wavefunctions and of the coeffi¬ 
cients of the linear combination is out of the scope of 
this work. Instead, following Ref. we model the exci¬ 
tonic wave function using the DFT wavefunctions of the 
lowest conduction band for the electron (resp. highest 
valence band for the hole) and take £k=k i V’kc'0ku as our 
excitonic wavefunctions, with the sum limited to a few 
representative points k, . We then approximate the exci¬ 
tonic lifetime by the sum of the electron and hole lifetime, 



FIG. 5: GaAs. Broadening (meV) of the critical point E\ as a 
function of temperature. Points: experimental data of Ref. (7; 
obtained by fit with the Fano-type excitonic line^. Squares: 
experimental data of Ref. 0 , fit with the two-dimensional crit¬ 
ical point model^. Solid black line: this work, Wannier inter¬ 
polation method extended to polar-optical coupling. Dashed- 
dotted line: theoretical calculation from Ref. 0 with empir¬ 
ical pseudopotential method. Both theoretical broadenings 
are due to electron-phonon coupling only and are calculated 
as the sum of the broadenings of the lowest conduction and 
highest valence bands, averaged over the four points along 
the A direction as explained in the text. The inset repre¬ 
sents the same experimental and theoretical results of Ref. 0 
as the main figure, with the same notations. Grey solid line 
on the inset figure: this work, Wannier interpolation method 
without extension to polar-optical coupling. 



FIG. 6: GaAs. Broadening (meV) of the critical point 
E 2 as a function of temperature. Squares: experimen¬ 
tal data of Ref. 0 , fit with two-dimensional critical point 
model^. Dashed-dotted line: theoretical calculation from 
Ref. 0 with empirical pseudopotential method. Solid black 
line: this work, Wannier interpolation method extended to 
polar-optical coupling, sum of the broadenings of the lowest 
conduction and highest valence bands, averaged over three 
points (see text). 
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r kc + r kl) . For the critical point E\, the width r kc (resp. 
r k „) of the electronic (resp. hole) level are obtained via 
eq. TO with four points k, equal to L,jL,^L and jL. as 
done in Ref. @. In the case of E 2 , only one representative 
point k = 2 - 7 r/a(|, \) was used in Ref. The region 

in the k space where valence and conduction bands are 
parallel and which contributes to the E 2 point was de¬ 
scribed in Ref. [43]. In this work, we decided to take into 
account three points k,: (§, 3 ), (1, |, |) and U which 

belong to the region which contributes to E 2 ^. 

1. Broadening of the E\ critical point 

The resulting broadening T kc + r k „, is reported as a 
function of temperature for the critical point E\ (Fig.0. 
As one can see, broadenings of the critical point Ei cal¬ 
culated in this work are in satisfactory agreement with 
the experimental results of Ref. In the case of our 
calculation, the agreement is best with the experimental 
data obtained by fit with the Fano-type excitonic shape, 
whereas the previous theoretical result of Ref. @ privi¬ 
leged the fit with two-dimensional critical point model^. 
The question to discriminate between the two methods of 
fit of the experimental data is, however, out of the scope 
of present work. Indeed, here we only demonstrate that 
the calculated widths due to electron-phonon scattering 
allow to correctly estimate the magnitude of the broad¬ 
enings of critical points, within the experimental error 
bar. 

In the inset of Fig. [5j we show the result of the calcula¬ 
tion of the same broadening of Ei , but with the standard 
interpolation method. As one can see, the broadening of 
E\ is slightly lower if the Frohlich coupling is omitted, 
however, the overall result is very similar, confirming 
the statement of the previous section that the scatter¬ 
ing of the electrons/holes away from valley minima is 
dominated not by the Frohlich, but by the intervalley 
scattering. 

2. Broadening of the E 2 critical point 

In Fig. [HI we have estimated the broadening of the crit¬ 
ical point E 2 in GaAs. In the case of E 2 , our results (solid 
line) are very different from the theoretical result ob¬ 
tained with the empirical pseudopotential method^, and 
yield a much better agreement with experiment, show¬ 
ing that the experimentally measured broadening of the 
E 2 point can be attributed to the electron-phonon scat¬ 
tering. It must be noted down, however, first that only 
experimental results extracted with the fit with the two- 
dimensional critical point model are available in this case. 
Second, the experimental behaviour of the broadening 


beyond 500 K differs from the one predicted by the cal¬ 
culation. 

In dotted line, we show the broadening of the critical 
point E 2 estimated in our work with only one point k = 
27r/a(|, |, |) as was done in Ref. [jl As one can see, the 
broadening reported with dotted line is similar to the one 
obtained with three representative points, showing that 
the result does not depend crucially on the method of 
averaging and that it remains widely different from the 
one obtained with empirical pseudopotential. 


IV. CONCLUSION 

In conclusion, we have presented the description of 
the extension to polar-optical coupling of the method 
which allows to interpolate the electron-phonon matrix 
elements in the space of maximally-localized Wannier 
functions. The extended method is based on Vogl’s model 
of the Frohlich component of the electron-phonon cou¬ 
pling, and allows to interpolate the electron-phonon ma¬ 
trix elements calculated within DFPT in polar semicon¬ 
ductors with excellent precision. We have applied the 
extended method of interpolation in the case of GaAs, 
and calculated the widths of the electronic levels due 
to the electron-phonon coupling for highest conduction 
and lowest valence band. We have demonstrated that 
the obtained widths of the electronic levels can be used 
to estimate the relaxation times of hot electrons and 
the broadenings of the critical points due to electron- 
phonon scattering, in good agreement with various ex¬ 
periments. Finally, we have shown that, although the 
Frdhlich interaction is the dominant scattering process 
for electrons/holes close to the valley minima, in agree¬ 
ment with low-field transport results, at higher energies, 
the intervalley scattering is expected to dominate the re¬ 
laxation dynamics of hot electrons or holes. 
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